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Abstract. Inverse scattering theory is extended to one-dimensional Schrodinger 
problems with near-boundary singularities of the form v(z —>■ 0) ~ —z~ 2 /A + 
v—\z~ 1 . Trace formulae relating the boundary value vq of the nonsingular part 
of the potential to spectral data are derived. Their potential is illustrated by 
applying them to a number of Schrodinger problems with singular potentials. 


1. Introduction 

Inverse scattering theory, developed more than 60 years ago by Gel’fand, Krein, 
Levitan, Marchenko, and others m la El SB El, enables one to reconstruct potentials 
of Sturm-Liouville or Schrodinger problems from scattering data provided they satisfy 
certain mathematical conditions In the case of Sturm-Liouville problems on the 
half-line, the potential v(z) is typically required to be continuous, to decay sufficiently 
fast at infinity so that it is absolutely integrable, and to remain finite on approaching 
the boundary z = 0. In a recent paper [7] (hereafter referred to as I), we investigated 
the exact solution of the 0(n) </> 4 theory on a three-dimensional film bounded in one 
direction by a pair of planar free surfaces in the many-component limit n —> oo. As 
has been recognised a long time ago in studies of the semi-infinite case where one 
surface plane is located at z = 0 while the other is at z = oo, the exact n = oo 
solution of this model at the bulk critical temperature T c leads to a one-dimensional 
Schrodinger equation with a singular potential v(z) = —z~ 2 /4, [8][9]. For temperatures 
above and below T c , one must deal with (rescaled) potentials that minimise a free- 
energy functional. The solutions of the corresponding equations for these potentials 
involve the eigenfunctions and eigenvalues of the Sturm-Liouville problem with 
the same potential. In other words, the potentials are self-consistent solutions of 
these equations. Following a standard practice in the physics community, we will 
occasionally indicate this fact by referring to them as “self-consistent” potentials. 

Unfortunately, the self-consistent potential that gives the exact solution of 0(n) 
<^> 4 model in the limit n —> oo is known only in numerical rather than in closed 
analytical form, barring the semi-infinite case at T c HDuniiiaiisi. From consistency 
requirements with the behaviour at T c and results obtained via boundary-operator 
expansions mmmmm one can conclude that they must vary as 

v(z ->■ 0) = + v 0 + 0{z) 

f For background on inverse scattering theory and an extended list of references, see [6]. 


(i.i) 
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on approaching the 2 = 0 surface plane and have a similar singular near-boundary 
behaviour at the second plane. 

The aim of this paper is to explore how standard tools of inverse scattering theory 
such as the Povzner-Levitan representation and the Gel’fand-Levitan equation [B] 
can be extended to one-dimensional Schrodinger problems on the half-line and finite 
intervals with potentials that exhibit the singular behaviour ED at the boundary. 
In addition, we will derive trace formulae for Sturm-Liouville problems with such 
potentials whose combination with inverse scattering tools will enable us to express 
the potential coefficient Vo through scattering data. In I, we applied some of the results 
obtained here to the analysis of the exact n —> oo solution of the above-mentioned 0{n ) 
</ 4 model, using them to determine a number of quantities exactly. Specifically, we 
there used trace formulae to be derived below and their corollaries to determine the 
exact value of the expansion coefficient i>o in equation (11.11) . Although trace formulae 
of a similar kind have been discussed in the literature both for nonsingular potentials 
as well as certain types of singular potentials, no trace formula for the case of potentials 
with the particular singular near-boundary behaviour ED was available. 

To put things in perspective and clarify the basis of knowledge from which we can 
start, it will be helpful to briefly recall the relevant literature. Singular potential terms 
of the form nz~ 2 with /r = l(l + 1), l = 0,1,... naturally arise through the centrifugal 
contribution to the effective potential of the reduced radial Schrodinger equations 
of radially symmetric three-dimensional Schrodinger problems. The corresponding 
inverse scattering problems with nonnegative integer values of the angular momentum 
quantum number l have been extensively studied for more than 6 decades, see [6] 
and its references. More recently, Sturm-Liouville problems pertaining to radial 
Schrodinger equations with non-integer values of l have also attracted considerable 
interest mmmmm- Freiling and Yurko [18] investigated inverse spectral 
problem for the case of non-integer l > —1/2. They described how the potential can 
be recovered from the spectral data, and proved a uniqueness theorem. Potentials with 
boundary singularities of the form specified in equation (11.11) (corresponding to the 
radial Schrodinger equation with l = —1/2), have previously been considered by two 
other groups of authors. Zhornitskaya and Serov m established a uniqueness theorem 
for the corresponding inverse spectral problem with vanishing potential coefficient 
V-i = 0. Kostenko, Sakhnovich and Teschl m succeeded to do this for the more 
general case V-i ^ 0. Unfortunately, neither of these two papers clarified the question 
of how the considered singular potentials can be recovered from the spectral data. 

Even if the potential does not become singular at the boundaries, the process of 
recovering it from scattering or spectral data is not easy because it usually requires 
the solution of integral equations [6]. If the potential v{z) is analytic near the origin 
2 = 0, one can try to solve the simpler problem of determining one or several leading 
contributions of the potential’s Taylor expansion about 2 = 0 from the scattering or 
spectral data. 

The first “trace formula” of this kind was obtained by Gel’fand and Levitan [2], 
who expressed the boundary value of a regular potential of the finite-interval Sturm- 
Liouville problem in terms of the corresponding eigenvalues. This seminal paper has 
triggered a lot of interest in such trace formulae owing to their rich applications 
in different fields of mathematics and physics, among them integrable nonlinear 
Hamiltonian systems. The Gel’fand-Levitan trace formula [2] was generalised in a 
variety of ways by Dikii [22] , Newton [23] , Faddeev and Buslaev [24] [25] , and others. 
An extensive list of references can be found in [6] EH [27], [28] - The trace formula 
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that comes closest to the one we need and shall derive below is the following one due 
to Newton [23] and Faddeev [MJ , which relates the expansion coefficient u 0 of the 
Laurent series 

vi(z —> 0) = ——-—- + vo + 0 ( 2 ) (1-2) 

z A 

of the effective potential vi of the radial Schrodinger equation on the half-line 
0 < z < 00 for angular momentum l to the scattering phase rji(k) and the discrete 
energy levels e„. It reads (see equation (VII.4.2) of reference [B]) 

8 f°° d 4 

= ~W+7j* Jo ( L3 > 

However, there are two reasons why this formula, equation ed, cannot be applied 
to the case of a potential with the boundary singularity (ED: The first is that the 
right-hand side of (11.31) does not exist for l = —1/2; the second is that the first-order 
pole term v-i/z of the potential ED is absent in ED- 

The remainder of this paper is organised as follows. In the next section we begin 
by defining more precisely the type of Schrodinger problems we will be concerned with. 
Both Schrodinger problems on the half-line (0, 00 ) and on a finite interval (0, N) will 
be considered. We then recall some necessary background of inverse scattering theory 
for nonsingular potentials, discuss its generalisation to potentials with near-boundary 
singularities of the form specified in equation ED, and explain the modifications of 
the required equations of inverse scattering theory implied by these singularities. In 
Sec. [3l we turn to the derivation of trace formulae for Schrodinger problems with 
singular potentials on finite intervals [0,1V] and the half-line [0,oo). Both the trace 
formula and its corollary we exploited in I are proven. In Sec. [4] we demonstrate 
the potential of these trace formulae by applying them to a number of illustrative 
examples. Section [5] contains concluding remarks. 

2. Schrodinger problems for potentials with singularities at boundaries 

We will be concerned with the differential equation (Schrodinger equation) 

/ H v ip(z) = e ip(z) (2.1) 

where TL V denotes the operator 

H v = -<9? + v(z). (2.2) 

Here 2 ranges either over the half-line (0,oo) or the finite interval (0,1V). Note that 
both z and N are dimensionless variables: they correspond to the distance from 2 = 0 
and the length of the interval measured in units of an appropriate reference length 
scale, respectively!! 

In the half-line case we assume the potential to have a Laurent series about 2 = 0 
of the form 

v(z) = v ( ~ sg \z) + u(z) (2.3) 

§ In I, a number of distinct reference lengths scales were encountered, namely, the lattice constant a 
of the discrete model which we used as a starting point there, the correlation and Josephson coherence 
lengths l/|m| (depending on whether the temperature variable m was ^ 0 ), and the film thickness 
L = Na. In applying the results of the present paper to I, the proper identification of the variables 2 
and N must be made. For example, depending on whether \m\ or L is scaled to unity, 2 : corresponds 
to the variables denoted z and 3 , respectively, in I. 
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with a singular part 

v (sg) (z) = + v-i\ (2.4) 

and a nonsingular contribution (Taylor series part) with the series expansion 

OO 

u ( z ) = '52 v j z3 - ( 2 -5) 

j=o 

For the sake of simplicity, we take u(z) to be analytic for all z > 0. Although 
the self-consistent potentials that yield the exact large-n solution studied in I are 
expected to involve non-analytic contributions of order z 3 In z, this is not a severe 
restriction and could be relaxed because we shall essentially rely only on the existence 
of u( 0) = vo and its first derivative u'(0) = V\. We also assume that u(oo) = 0 and 
that v(z) vanishes sufficiently fast for z —> oo so that 

pOO 

/ dz \ v{z)\ < oo for all zq > 0. (2.6) 

J z 0 

We shall also consider the analogous problems on the finite interval / = [0, iV]. 
When dealing with this finite-interval case, we shall assume the potential to have the 
symmetry property 

v(z) = v(N — z) (2.7) 

and to be decomposable as in equation (12.31) . Note that because of the symmetry (12.71) . 
u(z) contains singular boundary terms of the form v ( sg ) (N — z) unless it is explicitly 
stated that no singular part is included in the potential. Depending on the choices of 
I = [0, oo) and I = [0, A7] and the presence or absence of a singular potential term, 
we will refer to these cases as half-line or finite-interval problem with singular or non¬ 
singular potentials, respectively. To fully define these problems, we must also specify 
the boundary conditions. 

Consider first the half-line case with non-singular potential. It will be appropriate 
and sufficient for our purposes to choose Dirichlet boundary conditions ^(0) =0 and 
assume the absence of bound states. Extensions to Neumann and Robin boundary 
conditions are straightforward and are described, for example, in |4j and [Bi. Likewise, 
generalisations to cases in which bound states appear can be found in these review 
articles and their references. 

Let e = k 2 and <p(z, k) and <p(z, k) with £ = k 2 denote solutions of this Dirichlet 
problem on the half-line with two non-singular potentials v(z) = u(z ) and v(z) = u(z), 
respectively, which are normalised such that 

d z cp(z,k) | z =o = 1, d z ip(z,k) | z=0 = 1. (2.8) 

These regular solutions are related by a unitary transformation, which entails the 
so-called Povzner-Levitan (PL) representation [H [5} |6j 

tp(z, k) = ip(z, k) + ( dz' P(z, z') ip(z',k). (2-9) 

Jo 

Here the PL kernel P depends on the potentials v and v, but not on the spectral 
parameter e. From rigorous mathematical work (see 0 ), it is known that the kernel 
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satisfies the conditions [jj 

p ( z > z ) = \J ds[w(s)-w(s)], (2.10) 

P(z, 0) = 0, 

and the partial differential equation 

(5 2 - d 2 )P(z, s ) = [v(z) - 0(a)] P(z, s), 0 < s < z. (2.11) 

A case of particular interest is the choice v = 0. In this case the regular solution 
becomes 

sin(fcz) 


ip{z,k) = 


( 2 . 12 ) 


As is known from the classical work of Gel’fand and Levitan pH (see also 0) , one 
can reconstruct from the scattering data the function 


K(z,z') = j y(z,k)ip(z',k)[dp(k 2 )-dp(k 2 )\, 


as 


c\k k 2 

K(z,z') = 2 - ip(z,k)<p(z',k) 

Jo n 


U F(k)\ : 


- 1 


(2.13) 


(2.14) 


where F(k ) is the Jost function, while dp(e)/de and dp(e)/de denote the densities of 
eigenvalues e = k 2 of the Hamiltonians H v and Hy, respectively. The function K(z, z') 
can then be used as input for a Fredholm integral equation, the Gel’fand-Levitan (GL) 
equation 

K(z,z') + P(z,z') + [ dz" P(z, z") K(z 1 , z") = 0, (2.15) 

Jo 

from whose solution for the PL kernel P{z, z') the potential v{z) can be recovered via 


v(z) = v(z ) + 2 — P(z, z). 
dz 


(2.16) 


The above equations carry over to the Dirichlet problem on the interval [0, Af] with 
non-singular potential, with obvious adjustments. The spectrum becomes discrete. 
Let £j, = k 2 and e v = fc 2 denote the eigenvalues one obtains for the potentials v(z) 
and v(z), respectively, and ip(z,k u ) and Jp(z,k„) be the associated regular solutions. 
Instead of equation (12.141) . we then have 


K(z,z') = 


V=\ 


Jp(z,K)Jp{z',k v ) ip(z, k^ifiiz', k v ) 


(2.17) 


where a v and k v are the squares of the L2([0, N]) norms of ip(z, k v ) and ip(z, k v ), e.g., 
r N 


dz \ip(z,k„)Y. 


(2.18) 


Both the GL equation (12.151) and equation (12.161) remain valid. 

We next turn to the cases with singular potentials. We wish to apply the above 
equations (12.101) . (12.111) and (12.141) (12.161) to problems with potentials v(z) of the form 
specified by equations (12.31) (12.51) and their analogs v(z) = v^ se \z) with a vanishing 


|| The condition (12.10D holds when either Dirichlet or Neumann boundary conditions are applied. It 
gets modified in the case of Robin boundary conditions [4] . 
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nonsingular part. Whether these equations carry over to potentials with singular parts 
is a priori not clear. The original (non-rigorous) derivation of the GL equation in [I 
at first glance may seem to be generalisable in a straightforward fashion to cases with 
singular potentials. However, on second thought one realises that subtle issues such 
as the choice of boundary conditions and the interrelated ones of the existence of self- 
adjoint extensions of the Hamiltonian H v (|2.2D and of the completeness of the set of 
eigenfunctions arise for potentials v(z) involving singularities of the form (12.411 . 

These issues have been investigated in some detail for the critical-point 
Hamiltonian 

Uc = ~ h (2 ' 19) 

over the interval [0,1] in a recent paper [2!5|. But this work does not fully cover 
the cases we are concerned with. First of all, no additional singular term V-\/z was 
included. Second, no contribution that becomes singular at the plane z = 1 was 
included. Nevertheless, the paper provides helpful guidance. Let us briefly summarise 
some of its results relating to our investigations. 

The authors (KLP) of 129] show that the boundary conditions which square- 
integrable functions if belonging to the maximal domain T’maxl’Hc) = {if £ L 2 ([ 0,1]) | 
H c if £ L 2 {[ 0,1]) must satisfy at z = 1 can be parametrised through an angle j3 2 £ [0,7r) 
such that they can be written as 

if'(l) cos/32 + V’(l) sin ^2 = 0. (2.20) 


The choices /3 2 = 11/2 and /?2 = 0 correspond to Dirichlet and Neumann boundary 
conditions, respectively. 

Owing to the singular behaviour of T-L c at z = 0, the boundary conditions at this 
plane must be defined via a limiting procedure. KLP also show that any function 
£ Hmax(Hc) can be decomposed as 

if{z) = [ci(if) + c 2 {if)hiz\^fz + if{z), ( 2 . 21 ) 

where Cj(if), j = 1,2, are constants and if is a continuously differentiable function 
with Hcif £ L 2 ([0,1]), if(z) = 0(z 3 / 2 ), and if'(z) = 0(z 1 / 2 ). This leads them to 
consider self-adjoint realisations of H c on subspaces 


^ 01,02 


£ AnaxCHc) 


cos/3i Cl (VO + sin (3ic 2 (if) = 0,1 
cos/ 3 2 V>(1) + sin/3 2 V’(l) = 0 /' 


( 2 . 22 ) 


which are maximal in the sense that 


{if £ ^max {U c )\{'H c (f\if) = {(f\H c if) for all cf £ Vp 1} p 2 ] = T>p lt p 2 . 

According to KLP, given any f3 2 £ ]0,7r), one can choose anyone of these maximal 
subspaces 'Dp 1 ,p 2 as domain of H c and obtain thereby a self-adjoint realization, where 
the choice fS\ = 7t/ 2 corresponds to what is known as Friedrichs realization. 

For the special case of the Friedrichs realization, i.e. /3i = n/2 and hence 
C 2 (VO = 0, KLP found that the trace of the resolvent Tr(i? — He) -1 , the heat 
kernel Tr(e _tWc ), and the zeta function C(s,H c ) = Tr("H“ s ) exhibit a “usual” type 
of behaviour of the following kind: If \E\ —> 00 , with E belonging to any sector (of 
given opening angle) that does not intersect the positive real axis, then the first two 
quantities have series expansions in powers of (— E) k / 2 and ^( fc — 3 )/ 2 , k = 1 , 2 ,..., 00 , 
respectively. Furthermore, the zeta function C(s, H c ) can be continued from Re s > 1/2 
to a meromorphic function on the complex plane C with poles at s = 3/2 — fc, 
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k = 1,2,..., oo. However, for all other self-adjoint extensions (fli € [0,7r) with 
/?! y^ 7 t/ 2), logarithmic anomalies appear in the corresponding expansion of these 
quantities and t^(s,R c ). Thus the latter zeta function have a logarithmic branch point 
at s = 0 for such self-adjoint realisations. 

Let us now see how these findings relate to our Schrodinger problem with singular 
potentials on the interval [0, AT]. We are interested in the case where boundary 
conditions of the form specified in equation (12.2111 with C 2 = 0 hold at z = 0 and 
z = N. Owing to the symmetry of the potential with respect to reflections about 
the mid-plane z = N/2, the operator 7~L V commutes with the corresponding reflection 
operator R mi d ■ ip(z) i-» ijj(N-z). We can choose eigenfunctions <p v (z) = ip(z, k v ) that 
are even with respect to these reflections for odd values v = 1,3,..., oo, but odd with 
respect to R m id for even values v = 2,4,.. ., oo, by imposing the boundary conditions 

<p(z,k v ) = yfz[l + 0(2)], (2.23a) 

-2—^0 + 

ip{z,K) = (-l^v/F^K [1 + 0(iV - z)]. (2.236) 


The even eigenfunctions then satisfy Dirichlet boundary conditions at the mid-plane 
z = N/2, ip v {N/ 2 ) = 0, v = 1 , 3 , 5, ...; the odd ones Neumann boundary conditions 
ip' v (N/2) = 0, v = 2,4,.... The Hamiltonian R v becomes block-diagonal when 
represented in this basis. Hence, for v(z) = v^ ss \z) with v-\ = 0, we are back 
to the situation with /3\ = 7r/2 investigated by KLP, and we are dealing with the 
Friedrichs extension of the even and odd ‘blocks’ of R v . This means, in particular, 
that in the case of a potential which involves only the leading singular term specified 
in equation (12.41) (i.e., V-\ = 0 and no nonsingular contributions), the trace of the 
resolvent, the heat kernel, and the zeta function </(s,R v ) should display what KLP 
termed the “usual” behaviour and hence not involve logarithmic anomalies. 

However, it is clear that the addition of singular terms of the form v—\/z will 
destroy this usual behaviour and lead to the appearance of logarithmic anomalies. 
As we know from previous work mmm and can also be seen from the results 
discussed in Sec. VII of I, the surface free energy f s has a leading thermal singularity 
~ m 2 ln|m| whose proportionality coefficient is oc u_i , where it should be recalled 
that m oc (T — T c )/T c . Consequently, logarithmic anomalies do occur when v-± y^ 0 
in quantities such as the surface free energy and the surface excess energy. Owing to 
the close relationship of the above-mentioned three quantities investigated by KLP 
with quantities such as the energy and free energy of our model, similar logarithmic 
anomalies are to be expected in the former when u_i 0. Although a generalisation 
of KLP’s work to singular potentials (12.41) with v_i 0 appears possible, no such 
attempt will be made here since this is beyond the scope of the present paper. 

With regard to our subsequent considerations let us simply stress that the self- 
adjoint realisations of R v and we are concerned with are specified by the 

boundary conditions (12. 23 all and (|2.23fel) and their analogs for ip(z,k), respectively, 
where we assume that both sets of functions {g>(z, k^)}//L 1 and {Ip(z, k v )}//L 1 satisfy 
standard completeness relations, viz., 


E 


ip(z,k v )<p(z',k v ) 

k v 


S(z — z'). 


(2.24) 


We can now follow the steps taken by Gel’fand and Levitan in their original 
derivation [I] of the GL equations to convince ourselves that the PL equations hold 
also in our case with singular potentials. The starting point is the observation that 










Inverse scattering theory and trace formulae... 


the PL representation (12.91) for ip{z, k ) can be inverted to obtain a corresponding one 
for if(z,k), 

<p(z,k) = ip(z,k) + f d z' P(z,z')(p(z',k) (2.25) 

Jo 

with a different kernel P(z,z'). Here the regular solutions tp(z,k) and fi{z,k) for 
arbitrary k > 0 satisfy the boundary conditions 

= Vz[l + 0{z)}, 

z— >-0+ 

$(z,k) z = 0+ y/z[l + O(z)\. (2.26) 

Equation (12.251) tells us that ip(z,k) depends for given z £ (0,1) on all ip(z',k) with 

0 < z' < z. 

Following reference pQ, we can now choose z and z' < z in the completeness 
relation (12.241) . Then the right-hand side vanishes. Let us introduce the scalar 
product (X,Y) = X V Y V between infinite-dimensional vectors X = {X u ) 

and Y = ( Y v ) with real-valued components, using the coefficients 1/xy defined in 
equation (12.181) as (fixed) metric coefficients. The vanishing of the left-hand side of 
equation (12.241) then means that the vectors <l>(z) = (y>(z, k v )) and $>(z') = {tp{z', k„)) 
are orthogonal to each other, ($(z), $(*')) = 0. Since we know from equation (I2.25[) 
that tp(z', k u ) depends only on ip(z", k v ) with z" < z', the vector (tp(z', k v )) must also 
be orthogonal to $( 2 ). Hence we have 


E 


ip(z,kv)fi(z',kv) 


= 0. 


Substitution of the PL representation (12.91) then gives 


(2.27) 


E 


<p{z, kv) ip(z', k v ) 
Kv 


r dz"p( z ,z") 
Jo 


”'E 

V=1 


&{z", kv) ip(z\ kv) 


= 0. 


(2.28) 


Both contributions to the left-hand side diverge in the limit z' —»• z. We can get 
rid of this divergence by subtracting the completeness relation (12.241) for the set of 
functions {(p(z',k v )}. This yields the GL equation (12.151) with the kernel K given by 

<233- 

In order to show that the differential equation (12.111) for P(z , s) remains valid, 
we proceed along the lines of Rundell [3] in his derivation of his equation (6.6): We 
subtract the Schrodinger equations , H v ip(z) = £<p(z) and T-L{,ip(z) = elp{z) from each 
other and use the PL representation (12.91) to express y"{z) and v(z) tp(z) as 

= JZ2 ${ z ) + i p ( z 1 z ) z )] + -P (1,0) & z ) <f( z ) 

+ [ ds p( 2,0 \z, s) <p(s) (2.29) 

Jo 

and 


v(z) <p(z) = v(z) ip(z) + [v{z) — u(z)] ip(z) + 


ds v(z) P(z, s) <p(s), 


(2.30) 
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where P^’ k \z, s ) = d J z d k P(z, s). We then integrate the term of e (<£—<£) = 
involving second derivatives twice by parts. We thus arrive at 
,d P(z,z) 


0 = 


dz 


- [v{z) - v(z)\ 


(fi(z) 


+ [P{z, s) d s fi(s) - P ( °’ 1} (z, s)fi(s)] s=0+ 

+ [ ds<p(s){p( 2 ’°\z,s)-P {0 ’ 2 \z,s) 

Jo 

- [v{z) -v{s)] P(z,s)}. (2.31) 

If we set e = k 2 , then <p(z) becomes the corresponding eigenfunction Ip v {z) = ip(z, k). 
These eigenfunctions should be complete for the chosen boundary conditions (12. 23 all 
and (12.2361) . In order that the term P(z,s)d s <p v (s)\ a -o + exists, we must have 
P(z,s) = 0(s 1 / 2 ) as s —► 0. We therefore require that P satisfies the boundary 
condition 


P(2,s) = V^[Q(*) + 0(«)]. (2.32) 

This ensures that the contribution [.. .] s= o+ in the second line of equation (12.311) 
vanishes. Hence, both the coefficient of <p{z) in the first line and the one of ip(s) in the 
integrand must vanish, so that the condition (12.161) and the differential equation (12.111) 
follow. 

Unfortunately, finding an analytic solution to the GL integral equation (12.151) in 
the full square [0, N] x [0, N] is not normally possible except in a few special cases. Even 
in the semi-infinite case where K (z, z') can be expressed in terms of the eigenfunctions 
tp(z,k) and the Jost function according to equation (12.141) . analytic solutions do not 
appear to be feasible in general. 

If one were able to find analytic solutions for the two choices of potentials 
v(z) = v^ sg \z) + u{z) and v(z) = v ( ' sg \z) with u(z) of the form specified in 
equation m in the finite-interval or half-space case, then u(z) could be reconstructed 
via equation (12.161) as 

u(z) = 2— P(z, z). (2.33) 

dz 

Rather than aiming at a full analytic solution, wc will content ourselves here with 
a somewhat less ambitious agenda: using appropriate series expansion ansatze, we will 
determine the PL kernel P(z, s ) for small values of its arguments z and s. We first 
consider the simpler case in which the contribution V-i/z to v^(z) is absent and the 
Taylor expansion (12.51) of u involves only powers of z 2 . 


2.1. The case of v-\ = 0 and a nonsingular potential of the form u = u{z 2 ) 

We set V-i = 0 in (12.41) . together with all coefficients V 21-1 = 0 with l = 1,2,... ,00 
in the Taylor expansion m of u(z), and use the ansatz 

OO 

P(z,s) = y/zs ^ P hk Z 2 ’s 2k . (2.34) 

j,k =0 

Inserting both series expansions into equation (12.161) and equating the series -expansion 
terms of orders z°, z 2 , and z 4 yields the conditions 
i>o = 2P 0 o ; 
u 2 = 6(Pi j0 + Pop), 

V 4 = 10(P2 ! o + Pi,i + Pop)- 


(2.35 a) 
(2.356) 
(2.35c) 
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Likewise, matching the expansion terms of orders 

y/zs , Z 2 y/Is , 

and s 2 \fzs in 

eauation (12.111) rives 





O 

1 

4Po,i — P 0 ,oW 

= 0 , 


(2.36 a) 

16P2,o - 

- 4Pi,i — Pi.ouo — P 0 , 0 V 2 

= 0 , 


(2.366) 

4Pi,i - 

16Pq,2 - Po,iVo 

= 0. 


(2.36 c) 


These six equations for the six unknowns with 0 < j + k < 2 can be solved 
in a straightforward fashion. The results are 


p _ W 

Pro — — 


P-i.o = 


2 ’ 

12 16’ 
O 

V2 


PO.1 = TT 


12 16’ 


P - V4 
P.o — — 


60 384 


(2.37a) 
(2.37 6) 
(2.37 c) 
(2.37 d) 
(2.37 e) 

60 ' 384 96 ' ( ' 2 ' 37 ^ ) 

The corresponding expansion of the kernel K(z, s) can then be determined from 
the GL equation (12.151) . One obtains 

K(z,s) = V^\- V -Z+( V ^- V ^\ (z 2 + s 2 ) 


u o 


Pn = ^- 
’ 15 192 


P - v 4 

Pi,2 — — 


V, 


0 


V 0 V 2 
“32“’ 
VqV2 
48 ’ 

V 0 V 2 


, u 0 "2 \ / 2 
2 " ^ (l6 _ 12/ ( - 

_ _ufj_ _ 14 N , 4 , 4 n 
v 96 384 60y j 

(V0V2 , 

l 24 96 15 


(2.38) 


2.2. The case of v -1 ^ 0 and a nonsingular potential u(z) 


We proceed to the general case of a potential whose Laurent expansion about z = 0 
has the form 


v(z) = 


1 

4~2 


v-i 

z 


00 


+!>*"• 

3=0 


(2.39) 


We allow all coefficients Vj, j = —1,0,..., 00 , in this expansion to be nonzero and 
include a contribution V-\/z in v(z) = v^ ss \z). In this case, the ansatz (12.341) does 
no longer work and must be generalised. We use instead one of the form 

OO 

P(z,s) = Y J Pi{s/z)z\ (2.40) 

3=1 

where 0 < s < z < 1. We insert this ansatz along with the Taylor series m d into 
equation (12.111) . make the change of variable s —T p = s/z and Laurent expand in z. 
Matching the terms of order z _1 and z° gives us the differential equations 

(1-p- 2 ) [Pi(p)+4p 2 P"(p)] =0 


(2.41) 
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(1 - p-^v.^p) + ( p 2 - 1) P"(p) - 2 pP'(p) + l (9 - p - 2 ) P 2 (p) = 0. (2.42) 

Two linearly independent solutions of equation (12.411) are ^fp and ^fpXwp. The 
second one cannot contribute because it would be incompatible with the boundary 
condition (12.321) for P(z,s). The coefficient of the first is fixed by equation (12.331) . It 
follows that 

Pi(p) = ^VP- (2-43) 

We can now insert this result into equation (12.421) . The solution of this differential 
equation that is consistent with the boundary condition (12.321) and equation (12.331) is 
given by 

P*(p) = -4u 0 u_!) 2 Fi(- 1 / 2 , - 1 / 2 ; 1 ; p 2 ) Jp + V ° V ~ l( ^ + p) (2.44) 

where 2 F\{a, /?; 7 ; z) is a hypergeometric function. The homogeneous counterpart 
of the differential equation (12.421) has a second linearly independent solution 
G 2 ’2 (p 2 1 1 / 41 / 4 ) given by a Meijer G-function. It does not contribute to the solution 
because its behaviour 


^2,0 

( - r 2,2 



7/4, 7/4\ 
l/4,l/4) 


= —Vp 

p—> o 7r 


2 + In 


4J 


O(p 5/2 ), 


(2.45) 


is incompatible with the boundary condition (12.321) . 

Upon substituting the series expansion (12.401) along with a similar ansatz 


K (z^s) = J2 K j(s/z)z J (2.46) 

j=i 

into the GL integral equation (12.151) . one can solve for the lowest-order term of I\\ in 
a straightforward fashion. The result 

Ki (P) = ~y/P (2-47) 


yields 


K(z, 


s) 





(2.48) 


Hence we have 

v 0 = -2 


d K(z, z) 


dz: 


2=0 


—2 lim z 1 K(z,; 
2^0 


(2.49) 


3. Trace formulae, their proofs, and consequences 

In this section we will present and prove a trace formulae for Schrodinger problems on 
finite intervals and the half-line whose potential involves singular terms of the form 
specified in Sec. [ 2 ] [see equations (12.31) and (12.51) ] on approaching the boundary. We 
start by recalling some background. 

Trace formulae (see, e.g., i) allow one to relate properties of the potential to 
spectral or scattering data. Perhaps the earliest formula of this kind is the following 
theorem due to Gel’fand and Levitan [2]. 
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Theorem 3.1 (Gel’fand-Levitan trace formula) Let e v he the eigenvalues of the 
operator —d 2 + v(z) on the interval [0,1] subject to Dirichlet boundary conditions, 
where v(z) £ C^fO, 1] satisfies the condition dzv(z) = 0. Then the trace formula 

1 OO 

4 Mo) + u(i)] = J2 ( v2?r2 _ £ ") 

V=1 

holds. 

In I a trace formula for singular potentials of the kind encountered in the exact 
n — y oo solution of the 0{n) cj) A model on a film was needed to relate the coefficient vq 
of the nonsingular part of the potential to its scattering data. Here we derive a more 
general trace formula from which the one exploited in I follows in a straightforward 
fashion: 


Theorem 3.2 (Finite-interval trace formula) Let Hi be the analog of the 
operator defined by equation 12.2\) for the potential 

( v^(z) for 0 < z < N/2, 

' (3.1) 


viz = 


.,(sg) 


(TV - z) for N/2 < z < N, 


where v^ sg \z) is the singular potential Let e v = k 2 and <p{z,k v ), v = 

1 ,2,..., oo, denote the eigenvalues and eigenfunctions of Hi, that satisfy the boundary 
conditions h2.23a\) and 12.23b)) . Furthermore, let ip(z, k v ) and e v = k 2 be 
the corresponding eigenfunctions and eigenvalues of the singular boundary problem 
specified bv equations I2.1\) - \2.7j , \2.23a\ ) and H2.23b\) . Then the coefficient vq satisfies 
the relation 

r N/2 


VQ 


v=l 


= 2 V ( --+ — 

\ if.. TS. I N 


d z \y{z) 


- 


>(*)]. 


(3.2) 


where h v and x v are the squares of the L 2 ([0, N])-norms ( f 2.18)) of fp[z,k v ) and 
<p(z,k v ), respectively. 

To prove this theorem we start from equation (12.491) and use equation (12.171) for 
the kernel K(z, z ) to obtain 


K(z,z) = z^2 

z—tO z —' 


1 


1 

k v 


+ S{z) + 0(z 2 ) 


with 


S (z) = '^2^ 1 f> 2 {z,k v ) - y 2 {z,k v ) 


V=1 


(3.3) 


(3.4) 


Since the small- z behaviour of S(z) is controlled by contributions with v 1 , we 
can replace the quantities fc„, k v , and >c v by their asymptotic large-^ forms. They are 
conveniently expressed in terms of the momenta 

a„ = ir(v -l/2)/N. (3.5) 

We assert that the above quantities behave asymptotically as 


k„ — a v H- 

a v 


7^2 + [ln(4a y TV) + 7 e] \ + 0(a v 3 ), 


(3.6) 
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1 1 


K + a v I 21V 2 ' N 


r N/2 
“7 / dz 


u(z) — v^ c \z) 


+ -^-[ln(4aylV) + te] ^ + 0(a 1/ 3 ), 


(3.7) 


and 

=- 1 - 0(a„ 2 ), (3.8) 

TTCLi/ 

where qg is the Euler-Mascheroni constant. 

To show this, note first that the regular solution f>{z, k) that satisfies the boundary 
condition (12.23al) at z = 0 becomes for k = k u the eigenfunction pertaining to the 
eigenvalue fc 2 . It is given by 

<p{z, k ) = e 17r / 4 —]L= i J—2\kz) = \fz e lkz [|(1 + iv-i/k); 1; —2ifez], (3.9) 

V2fc 2fc ’ 

where M Ktfl (z) denotes a Whittaker M-function [30l [3T1 [32] . When v-i = 0, the result 
simplifies to 

ip(z,k) = y/zJo(kz), V-i = 0, (3.10) 


where J v (u) is a Bessel function of the first kind. 

Second, the eigenvalues follow from the Neumann and Dirichlet boundary 
conditions the functions ip(z,k v ) fulfil at z = N/2 for odd and even values of u, 
respectively. 

The boundary conditions at z = N/2 yield for k v the equations 


0 

which simplify to 

0 


M' .t._i J-NiK), 

l ~2TT ’ u 

v = odd, 

M j»_i J-Nik v ) = 0, 

1 2k ’ u 

v = even, 


MNfeu/2) - Nk„ Ji(Nk v / 2), v = odd, 

Jo(Nk v /2), v = even, 


(3.11) 


(3.12) 


when V-\ = 0 . 

We first consider the solution of these equations in the simpler case V-i = 0. 
The asymptotic expansions of the roots of Jo can be found in equation (10.21.19) 
of [52] . From them, the result given in equation (13.61) with V-i set to zero follows 
at once for even v. To derive the corresponding i?_i = 0 result for odd is, one can 
use the familiar asymptotic expansions of the Bessel functions. Equating the leading 
contribution oc kj cos(Nk v /2 — 7 r/ 4 ) of the asymptotic series for the first line of 
equation (13.61) to zero yields = a u + 0(1 /a„) for odd is. Including the next-to- 
leading term oc k v of this series, one can expand in k v — a v to obtain for the 
difference the result 1/(2 N 2 a u ) + 0(l/a 3 ) in accordance with equation (13.61) . 

To do the analogous computation when v-\ ^ 0, we need the asymptotic large- 
k behaviours of ip(N/2,k) and its derivative d z ip(N/2,k). These can be determined 
by using the asymptotic expansion of the function M_; Kj o(— iz) for z —> oo given in 
equation (13.5.1) of [30] and then expanding the truncated series to first order in V-i. 
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The calculations are conveniently performed with the aid of Mathematic^ . One 
obtains 





(2k z)~' 1K T (1/2) 
I\l/2- rci) 


exp 


•7° • 7F 

i kz — l-1- 

4 


7 TK 
~2 


1 >2 / i 1 + i/c 

4 l 2kz (2 kz) 2 


4 exp 

nk 
x cos ^ kz 


TTV-i 

4 k 

7T 

_ 4 ” 


1 


--—- + 

16 k 2 z 2 4k 2 z 

1 u_i 
8 kz 2k L 


0(fc- 3 ) 


7t 2 z;i 


16fc 4 
7e + ln(8/cz) 


(3.13) 


+ 0 (fc" 7/2 ). (3.14) 

The boundary conditions at z = N/2 imply for v = 1,3, ...,oo that the 
cos{.. -} Z =N/ 2 must vanish, and for v = 2,4,..., oo one finds the condition [tanh{...}+ 
0 (fc _3 )] z= jv /2 = 0, where the ellipses in both cases represent the argument of the cosine 
in equation (13.141) . From these conditions, the result for k v given by equations (13.61) 
and m follows at once. 

To derive equation (ED, we can use first-order perturbation theory in u = v — v^ 
to determine the energy e v = k 2 . This gives 


K = \j°k 2j r (dp v \u\(p v )lk v + ... = k v + 4- (<p v \u\<p v )/k v + O (k u 2 ) (3.15) 

with 


f N/2 

{<p u \u\ip„) = 2 / u(z) <p 2 (z, k) dz. 


(3.16) 


Taking into account the familiar identity cos 2 x = [1 + cos(2x)]/2, we conclude from 
the asymptotic expansion (13.141) that 


ip (Z , kjy) - 


1 + cos [k„z + 7r/4 + 0(1 /k v )\ 


Ttkv 


1 + 0(fc“ 1 ) . (3.17) 


The integral cos[fcy^+0(l)]u(2;)dz vanishes as k u —> oo by the Riemann-Lebesgue 
lemma. Consequently, we have 


(Pv\u\pv) 


2 

7T k u 



u(z) dz + O (k v 2 ) 


(3.18) 


and likewise k v = N(k u n )~ 1 + 0(k~ 2 ). Equation (13.71) follows from these results 
in conjunction with equation (13.61) . Finally, equation (13.81) also follows because 
perturbative corrections of order u and higher are smaller by at least one factor of 
1 /a v . 


Wolfram Research, Computer code Mathematica, version 10. 
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Having established the above results, we can now exploit them to work out the 
small-z behaviour of S(z). An elementary calculation gives 

OO 

s(z ) = \a v (k v - k v )d av [<Pv(z, a v )] 2 + 0(a“ 1 ) + 0(z 2 ) 


U— 1 

= z lim 

z-> o 


[! + °K*)] d a v [${z, av)/y/z? 


v=l 


r i/2 

/ dz'u(z') + 0(z 2 ). 

Jo 


(3.19) 


Since [<£(z, a v )/^fz\ 2 is a function f(a v z), the required limit becomes that of a Riemann 
sum lim z _j.o7i f'(a v z) yielding the integral f//°dkf(k). Upon substituting 
v — v for u and exploiting the boundary condition (I2.23all at z = 0, we arrive 
at 

r-i/2 


S(z) = —z 




d z’ 


v(z')-v^(z') +0 (z 2 ), 


(3.20) 


whose insertion into equation (13.31) yields equation (13.21) and hence completes the proof 
of Theorem 13.21 

An immediate consequence of Theorem 13.21 is that the difference vo — vo of the 
potential coefficients vo and vo associated with two potentials v(z) = v^(z) + u(z) 
and v(z) = v^ c \z) + u(z) with the same singular part v^ c \z) can be written as 

r N /2 




dz [u(z) — u(z)], 


(3.21) 


where iq, is the analog of >c v . 

We next turn to the derivation of the trace formula for the half-line case used in 
I. It is not difficult to see that the following theorem can be inferred from Theorem l3.2l 
by taking the limit N —> oo. 

Theorem 3.3 (Trace formula for half-line case) Let v(z) and v(z) be two 
potentials on the half-line (0,oo) that vanish faster than z _1 as z —^ oo and behave as 

v(z —> 0) = y( ss \z) + vo + o(l), 

v(z — > 0) = i/ ss \z) + vq + o(l), (3.22) 

with the same singular part Let <p v (z), v = 1,..., nb, be the regular solutions to 

the Schrodinger equation 'H v (p v (z) = e v ip v (z) subject to the boundary conditions 

<p v (z) = y/z[ 1 + O(z)] (3.23) 

z—> 0 

that correspond to bound states (where nb may be zero). Denote by 

p OO 

x»= 'L'l(z) dz 


(3.24) 


the squares of the L 2 ([0, oo)) norms of these (real-valued) functions. Denote 
furthermore by A(k) = eand rj(k) the scattering amplitude and the scattering 
phase, respectively, which are defined through the large-z asymptotic behaviour of the 
regular solution 


<p{z,k) = 


A(k) 


sin[fcz + r/(k)\ + 0(1/z), 


(3.25) 
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in the absolutely continuous part k > 0 of the spectrum. Furthermore, let (p^, hb, 
and A{k) = e a ( k ' > be the analogous quantities pertaining to the potential v(z). Then 
the following relation holds 



pOO 

/ d k k 2 

e -2cr(fc) _ e -2 a(k) 

Jo 

■ 


+ E 

bound 
states v 


2 

Xv 


E 

bound 
states u 


2 


(3.26) 


To derive this theorem, note first that the integral on the right-hand side of 
equation (221) is 0(1 /N) and hence vanishes as N —> oo. Owing to the symmetry 
property urn the eigenfunctions for finite N come in pairs that are even or odd with 
respect to reflections about the mid-point N/2. For each even (odd) eigenfunction that 
becomes a scattering state in the limit N —> oo, there is an odd (even) eigenfunction 
approaching the same scattering state. For either one of these eigenfunctions of any 
such pair we can use the fact that the regular solutions can be approximated as 

<p{z, k„) w e a( ' kl, ' > sin[fcj,z + 7)(k v )\, 

p(z, k v ) w e CT( ^ ) sin [k v z + r){k v )\, (3.27) 

k v 


in the inner region < z < N — lof the film. The associated norm parameters behave 
as 


2 K 


(3.28) 


because the error resulting from the differences in the boundary regions is 0(1). Using 
this, one concludes that the resulting sum j? ... of those states v that turn into 
scattering states becomes an integral / 0 °°(dfc/7r) ..., which yields the integral on the 
right-hand side of equation (13.2611 . 

Next, consider the remaining states (if there are any), which become bound states 
when N —> oo. These also come in pairs of even and odd eigenfunctions that approach 
the same bound state on [0, oo). The factor of two is canceled by the corresponding 
factor of 2 in the norm parameters. Consequently, equation (13.261) follows, which 
completes the proof of Theorem 13.221 Hence we arrive at the following corollary 
(exploited in I). 


Corollary 3.1 Let v(z ) and v{z) be two potentials on the half-line (0,oo) with the 
following properties: 


(i) The Schrodinger equation TL v if{z) = stp(z) subject to the boundary condition 
if(z —> 0) = O (\fz) and its analog with v —> v have no bound-state solutions. 

(ii) Both potentials have the same singular behaviour at z = 0 specified in 

equation with identical coefficients V-i and V-\ though possibly different 

limiting values vq and vq °f their regular parts. 


Then the following relation holds between the difference of the latter coefficients and 
cr(k), the logarithm of the scattering amplitude: 


4 r 

vo ~ vo — ~ d kkr 
JO 


g-2o-(fc) _ g-2 a(k) 


(3.29) 
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To illustrate the potential of the above trace formulae and the Corollary 13. 11 we now 
discuss their application to several illustrative cases. Only potentials are considered 
whose leading near-boundary singularity — (4z 2 ) -1 for 2 —> 0 agrees with that of the 
self-consistent potential of the exact n —> oo solution of the 0(n) theory on a him 
discussed in m, m and I. 


4-1- Half-line case with v{z) = — (4sinh 2 2) 1 

We begin by discussing the application of Corollary 13.11 to the Schrodinger problem 
on the half-line with the potential 

1 


/%) = - 


4 sinli z 

From its Laurent expansion about z = 0, 


one finds 




„ (i) _ 1 
V ° 12 ' 


1 

4^ 


1 

12 


z 

60 


■ 0(2 


(4.1) 


(4.2) 


(4.3) 


In order to apply Corollary 13.ll we choose the auxiliary reference potential v(z) 


as 


u ( 1 ) (z) = - 


4z 2 ’ 


(4.4) 


so that the corresponding potential parameter vanishes. The same holds for the 
corresponding coefficients of the other auxiliary potentials v^\ j = 2,3,4, we shall 
consider below, i.e. 


-Xo) _ 


= 0, j = 1,2,3,4. 


(4.5) 


The regular solution satisfying the Schrodinger equation defined by equa¬ 
tions m, m, and 631 subject to the boundary condition (I2.23al) is 


<P { 1 \z,k) = \l- 


r(i/2 + ifc) C-i/ 2 (coth z ) 


= y/ tanh(z) cosh lk (z) 2 Fi ( 2lfc 4 +1 , 2lA ^ +,i ; 1; tanh 2 z) , 

(4.6) 

where Q0(C) an associated Legendre function of the second kind |30l 3T( [32]. The 
result given in the second line of equation (14.611 holds for z > 0 and k > 0 because 
Qv(() can be written as 

i 7T 1 / 2 r(/z + i/+l)(C 2 -l)' 1 / 2 


QS(C) = e""- 


■ 2 F 1 (^ + l,^±V + f;l/C 2 ) 


(4.7) 


2 v + 1 C+»+ 1 T(v + 3/2) 
when £ £ (1, oo). 

To determine the asymptotic large-z behaviour of </jW(z, k), we set t = 1 — £ = 
1 — cothz = 2e -2z + 0(e _4z ) and expand in t. We thus arrive at the limiting form 


V (1) {z,k) = 


F^\-k)e ikz - F^\k)< 
2i k 


—ikz 


(4.8) 
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with the Jost function 

= -i ■ 


(4.9) 


From the latter, the associated scattering amplitude A^(k) can be computed in 
a straightforward manner via the familiar relation A[k ) = y/F(k) F(—k) between the 
Jost function F(k) and the amplitude A(k). The result is 

A (1) (k) = e a(1){k) = [(2fc/7r) coth(fc7r)] 1/2 . (4.10) 


The regular solution for the auxiliary Schrodinger problem with the 
potential (14.171) problem is given by 

(p^\z,k) = y/z Jo{kz). (4-11) 

From its asymptotic behaviour 

(p^\z,k) = \ — sin(fcz + 7 t/ 4) + 0(l/z) (4-12) 

z->oo V 7TK 

we can read off the scattering data 

i (1) (z) = e* (fc) = y/2kj^, ?) (1) (fc) = (4.13) 

Upon inserting this result for a^(k) and its analog for a^\k) given in 
equation (14.1011 along with equation (14.511 into equation (13.291) . we recover indeed 
the value of vq stated in equation (14.31) : 

A r OC -i 

^o 1} = - J o dfcfc2 - tanh ( fc7r )] = J2' ( 4 - 14 ) 


4-2. Half-line case with v^{z) = —(4 z 2 ) 1 + 4(7t 2 z) 1 


As second example, we consider the Schrodinger problem on the half-line for the 
potential v^ 2 \z) = v(z; oo, 1) — 1, where v(z;L,m) represents the self-consistent 
potential of I. Recall that the exact analytical form of this potential is not known. 
However, its following properties established in I determine it uniquely: 

(i) The small-z behaviour of the potential v^ 2 \z) is given by equation EU with 

V-\ = 4 / 7T 2 . 

(ii) For z 1, the potential v^ 2 \z) vanishes exponentially, 

v^ 2 \z) = exp[— 2z + O(lnz)], z —>■ oo. (4-15) 


(iii) The Schrodinger operator on the half-line 0 < z < oo, defined by 

equation EH). does not have bound states. 

(iv) The scattering amplitude A^ (fc) determined by the large-z asymptotic behaviour 
(13.251) of the corresponding regular solution ip(z, k) reads 


A^ 2 \k) = y/k/ arctanA:. 


As auxiliary potential, we here choose 


v^ 2 \z) = 


1 4 

4z 2 tt 2 z ’ 


(4.16) 

(4.17) 


which fulfils equation ESI). As before, it is understood that the boundary 
condition (12.261) is imposed on the respective Schrodinger equations (12.11) . 
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According to I, the regular solution of the Schrodinger problem for the auxiliary 
potential is given by 

eiir/4 2 

<^ (2) (z, k ) = M_ iK}0 (-2ikz), k = 


V2k^- 1K ’" K ■' 

where M_j Ki o( z ) is a Whittaker-M function [301131 IS]. It behaves as 

4 k 2 


tp'- 2 '(z,k) = \fz 


2-S-O 


1 4 2 

1 H-+ 

7T 2 


0(2 3 


for small 2 . Its asymptotic large -2 behaviour is of the form 
<p( 2 \z,k) ~ ^ sin[fcz + rj^ (fc, z)\ 

z->oo k 

with the amplitude 

A^ 2 \k) = exp cr*' 2 '* (fc) = [l + e 4 /( 7rfc )] /7r 


and the logarithmic phase shift 


f] (2) {k,z) = A + argT ( ^ + 


2i 

7r 2 fc 


r 2 fc 


ln(2fc^). 


(4.18) 

(4.19) 

(4.20) 

(4.21) 

(4.22) 


The appearance of the 2 -dependent logarithmic term in fj(k, z ), the analog of the 
phase shift, is due to the slow decay ~ z of the potential for z —> oo and well-known 
from the case of scattering by a Coulomb potential. 

^ /q\ 

These results can be inserted into equation (13.2911 to compute Vq . Since the 
calculation of the required integral has been explained in I, we just quote the result 


,( 2 ) 


dfc k 2 


7r/fc 


arctan k 


1 + exp[4/(7rfc)] 


56C(3) 


- 1. 


(4.23) 


4-3. Half-line case with v^ ss \z) = —(4 z 2 ) 1 — 4(7t 2 2) 1 


We next consider the case of the Schrodinger problem on the half-line for the case of 
the self-consistent potential v( 3 \z) = v(z;oo,— 1). Again, the exact analytical form 
of this potential is not known, but the following properties known from I determine it 
unambiguously. 

(i) The small -2 behaviour of the potential v^ 3 \z) is given by equation CD with 
V-i = —4/7r 2 . 

(ii) In the limit 2 —> oo the potential ( z ) vanishes as 


t> (3) (2) 


1 

213’ 


2 —► OO. 


(4.24) 


(iii) The Schrodinger operator 3 ) on the half-line 0 < 2 < 00 , defined by 

equation (12.21) with the potential 1 ^(z), does not have bound states but a “half¬ 
bound state” [33] = <p(z,k = 0) with zero energy. The regular solution 

<po(z) approaches unity as 2 —> 00 , lim,-^ tpo(z) = 1. 

(iv) The scattering amplitude A/ 3 ) (fc) for this Schrodinger problem with the potential 

( 2 ) is given by 


A< 3 >(fc) 


k 

\/l + 7r|/c|/2 


(4.25) 
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As reference potential we here choose 


v {3 \z) = 


1 4 

4 Z 1 7 t 2 z 


(4.26) 


Thus, in the case of the Schrodinger problem with the potential v^(z), the 
spectrum of the Hamiltonian H v is continuous; no discrete (pure point) spectrum 
exists. However, the spectrum of the Hamiltonian T-Ly( 3 ) can be decomposed as 
spec('H S (3)) = { E „® [0,oo) into a pure point part with eigenvalues 


£0) = = _ 
V — 


tt 4 (t2- 1/2) 2 ’ 


12 = 1,2, 


(4.27) 


and a continuous part [0,oo) of improper eigenvalues k 2 indexed by k > 0. The 
eigenfunctions associated with the eigenvalues (14.271) are 


~( 3 )/ \ ^y/v-1/2 

‘Pv \ z ) = -o- 


= yfzeyi p 


M, 


4 z 

u - 1/2 ’°\ k tt 2 (u- 1/2) 

L 

L V -l 


8 z 


7T 2 (2:/- 1)J ~ 1/_i Ltt 2 (2^ — 1)J ’ ^ 4 ' 28 ^ 

where L„_i(z) are Laguerre polynomials. The regular solutions for k > 0 and k = 0 
are given by 


e W4 
\/2k 


<P {3) (z, k ) = M- wfi {-2\kz), k = 


(4.29) 


and 


(z, 0 ) = y/z J 0 ( 4 z 1 / 2 / 7 r), (4.30) 

respectively. Both c p^\z ) and (p^ 3 \z , fc) have been normalised such that the boundary 
condition (12.23al) is fulfilled. As shown in I, the large- 2 ; behaviour of the regular 
solutions (14.281) is of the form (14.201) with the scattering amplitude 

A^ 3 \k) = exp<j( 3 \k) = k [ 1 + e -4 /! 7 ^)] /n (4-31) 

and the logarithmic phase shift 

« <a)(t - -•) = 1 + " sr (5 - -§k) + A l,,(2kz) - (4 - 32) 

For the norm parameters k v we obtained in I 

7T^ 

* 1 /=-XT (2i2-i) 3 , 12 = i,2, ...,oo. (4.33) 

64 

The summation of the series and the computation of the integral one encounters when 
applying the trace formula (13.221) have been performed in I. We therefore again just 
quote the result 


i4 3) = - f d k 


nk 1 7 ik 

1 + exp[— 4 /( 7 rA;)] 2 


128 A 1 




7T 4 2^-2 (2v — l) 3 


56 C(3) 


so that Vq 3 ' 1 and v^ 2> + 1 have identical values according to equation (14.231) . 


(4.34) 
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To apply and check the finite-interval trace formula m , we now consider the 
Schrodinger problem on the unit interval (0,1) for the potential 

—.‘i 


v^\z) = -- 


4sin 2 (7rz) 
From its Laurent series about 0 = 0, 


we read off that 


i _2 „2 r 2 

,(4) M = _ L 7 f_J_ 

’ Az 2 12 60 


,( 4 ) 


+ 0( z4 )i 


U n = -- 


(4.35) 


(4.36) 


(4.37) 


0 12 ' 

It belongs to the class called “trigonometric Rosen-Morse potentials” in the literature 


The solution of our Schrodinger problem for this potential can be obtained with 
the aid of the solutions of the Schrodinger equation for the potential — U 0 / cosh 2 (ax) 
given in § 23 of [36]. Setting a = in, x = z — 1/2, and U 0 n 2 /4, imposing the 
boundary conditions, and taking into account that the eigenfunctions are again either 
even or odd with respect to the mid-point z = 1 /2, one can determine the eigenenergies 
E v and eigenfunctions tpv(z) in a straightforward fashion. One obtains 

E„ = [ 7T(*2 - 1/2)] 2 , 1/ = 1,2,..., 00, (4.38) 

and 

, x /sin(7 tz) r , 

Mz) = Y— ■^—^-P„-i[cos(nz)\, (4.39) 

where the Pj(x ) are Legendre polynomials. 

Recalling the familiar orthonormality property 

O' + V 2 ) / i da: P i 00 P 3' 00 = hr > ( 4 - 40 ) 

one sees that the squared norms of the eigenfunctions are given by 
r 1 i 

*v= \ ( PZ{z)dz=— - ——. (4.41) 

Jo n 2 (v — 1/2) 

As auxiliary potential we choose the potential v(z) defined in equation (13.11) with 
V—i = 0 and N = 1, namely 


for 0 < z < N/2, 

»(*»(,) = 4 -" , (4.42) 

( for N ! 2<z< N ' 

The eigenfunctions for this auxiliary Schrodinger problem satisfying the boundary 
conditions (I2.23al) and (12.2361) are even and odd under reflections about the mid-point 
0 = 1/2 when i/ = l,2,...,oois odd and even, respectively. They read 


y/zJo(kvz) 

(-l)" -1 yfl-z Jo[k„(l - z)} 


for 0 < 0 < i, 
for \ < 0 < 1, 


(4.43) 
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where the k v are to be determined from the Neumann or Dirichlet boundary conditions 
at z = 1/2, depending on whether v is odd or even. 

Let us use the standard notation j 0 ,m f° r the mth zero of the Bessel function 
Jo(k) and denote the mth zero of the function J 0 (fc/2) — k J\(k/2) as £ m , i.e., 

Jo(tm/2) - Im Jl(^ m /2) = 0. (4.44) 


Then we have 


kv 


Zjo,v/ 2 , v = 2,4, ...,oo, 
•^ 0 + 1 )/ 2 j v = 1, 3, • • ■, oo. 


For the squared L 2 ([ 0,1]) norms of the functions fiv{z) one finds 

k v = [J 0 (W2) 2 + Ji(fc,/2) 2 ]/4. 


(4.45) 


(4.46) 


The above results can now be inserted into the trace formula 13.21 to obtain 

vf = 2 E (£- ~ —) + 2 T d * [^ (4) (z) + Hz)] 

v—1 V Kv J Jo 

^- 1 = 0.822467033..., (4.47) 

where we evaluated the series Y^=\ numerically using Mathematica [?], verifying 
that the resulting number agrees to 9 digits with the exact value given in 
equation (14.371) . 



5. Concluding remarks 

In this paper we have extended standard tools of inverse scattering theory such as the 
Povzner-Levitan representation and the Gel’fand-Levitan equation to one-dimensional 
Schrodinger problems on the half-line and finite intervals whose potentials exhibit near¬ 
boundary singularities of the form m known from the self-consistent potential of 
the exact n —> oo solution of the 0(n) 4> A model. As we have seen, the boundary 
singularities of the potentials imply modifications of the analogs of these equations for 
nonsingular potentials. 

We have also derived new trace formulae for Schrodinger problems on the half-line 
and finite intervals for such singular potentials. As is borne out by these trace formulae, 
the boundary singularities entail important modifications of known trace formulae for 
nonsingular potentials (cf. [6]). The comparison of our trace formula (13.261) with the 
standard one stated in equation (11.31) reveals the differences. In the latter, the potential 
coefficient Vq is expressed in terms of the scattering phases r]i(k) and eigenenergies e v . 
By contrast, our trace formula (13.261) relates vq to the amplitude <r(k) and norm 
parameters n v . We checked these trace formula or the implied Corollary 13.11 here 
by applying them to four illustrative cases, for some of which exact results for the 
potential parameters Vq'\ j = 1,...,4, are available. In I we benefited from the 
trace formulae in that they enabled us to determine the potential parameter vo of the 
self-consistent potential. 

We expect both the inverse scattering techniques for singular potentials and the 
trace formulae described here to have useful applications to other interesting problems. 
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